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Abstract 

We generalize the well known formulation of the susceptibles, infected, 
susceptibles (SIS) spatial epidemics with creation and annihilation opera- 
tors to the reinfection model including recovered which can be reinfected, 
the SIRI model, using ladder operators constructed from the Gell-Mann 
matrices known in quantum chromodynamics. 



1 Introduction 

Stochastic epidemic models like the susceptibles, infected, susceptibles (SIS) 
spatial epidemics, also know as contact process, are generalized to include more 
realistic scenarios, e.g. including resistance and eventually milder reinfection, 
giving the SIRI stochastic model. An analysis in terms of pair approxima- 
tion has been recently given for the SIRI model to obtain phase transition 
lines |Stollenwerk, Martins and Pinto (2007)| . This improves the rough qual- 
itative analysis using mean field approximation, leading to the often studied 
simple ordinary differential equation systems. To obtain more insight into 
such stochastic epidemic models, strating from the master equation formulation 
van Kampen (1992)] a field theoretic formulation with creation and annihila- 



tion operators giving a Schrodinger like operator equation has been used, first in 



Doi (1976) , then subsequently elaborated [Grassberger and de la Torre (1979) 



Grassberger and Scheunert (1980) , Peliti (1985) and more recently used e.g. in 



Brunei, Oerding and van Wijland (2000) , Park and Park (2005)] and e.g. in 



de Qliveira (2006) 
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In the analysis of such stochastic spatial models Pauli matrices are used 
to construct spin 1/2 ladder operators which serve as creation and annihi- 
lation operators for a formulation of the master equation in vector notation 
Hinrichsen (2000)] for the SIS epidemic model. We generalize to the SIRI epi- 



demic model, finding ladder operators as creation and annihilation operators 
derived from the Gell-Mann matrices, normally used in quantum chromody- 
namics (see e.g. |Aguado (2008)| ). 

While mostly a subsequent renormlaization analysis based on path integrals 



is performed Brunei, Oerding and van Wijland (2000) , Park and Park (2005) 



so that classical methods from particle physics can be used [Zinn-Justin (1989) 



another way of series expansions based on a perturbation ansatz, as described in 
Dickman and Jensen (1991) , de Oliveira (2006)] is a promissing way to obtain 



e.g. phase separation lines in the SIRI model. Up to now only simple systems 
like the one dimensional contact process, i.e. the SIS epidemics, have been 
analysed in this way. 



2 The SIS epidemic model revisited 

The Pauli matrices give lowering and raising operators for spin 1/2 systems 
(named after Wolfgang Pauli, see Zinn-Justin (1989) |). The Pauli matrices are 

ff - := (i I) ' °"» := (i ~o) ' az:= { \ -l) (1) 

and from these the ladder operators Brunei, Oerding and van Wijland (2000)] 
are: the raising operator 



OA 



\{<? x + i<y y )=( ° I } V2) 



and the lowering operator 



Hence for the states 



|1>==(J) • *»==(!) , <«) 

we have creation operator c + = o+ giving 

c + |0) = |l) (5) 
and annihilation operator c = <7_ giving 

c|l> = |0> • (6) 
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Then the SIS epidemic with master equation 

N 

jp(h,...,i N ,t) = ^2 w hA - It (t) p(h,...,i- ii,...,i N ,t) 



dr 

i=i 



(7) 



N 



^2 w i-iijM p{h,-,Ii,-,lN,t) 



1=1 



for Ij 6 {0, 1} and transition rate 



N 



wi uX - h = -Vi } -h + a-il-Ii) , (8) 



and 



A' 



wi-k,U = P \J2 Wij -0--Ii) + <*'!* , (9) 
can be given in vector notation by 

±Mt))=L\9(t)) (10) 

for a Liouville operator L to be calculated from the master equation (Eq. ([7])) 
and with state vector \^{t)) defined by 

i*(«)> := £...5>(/i, ...,/*,*) K) 7l -(4) /Jv io) 

/ 1= /jv=0 

(11) 



{/} \i=l / 



|0) 



and vacuum state |0). 

Parameters of the model are (for birth-death process or SIS epidemic) (3 
birth or infection rate and a death or recovery rate. Here ( Jy) is the adjacency 
matrix containing for no connection and 1 for a connection between sites i 
and j, hence Jjj = Jji £ {0, 1} for i ^ j and Ju = 0. 

For lattice sites i we have creation and annihilation operators defined by 
c+|0) = |1) and a\l) = |0), and (c+) 2 |0) = and a\0) = 0, hence 

c 4 = h-\l-h) (12) 

c+|/,) = (1-70-ll-Ji) (13) 

and (c^) 2 |/j) = c||/j) = 0. We have anti-commutator rules on single lattice 
sites 

[ci,cf} + := acf + cfci = 1 (14) 
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and ordinary commutators for different lattice sites i ^ j 



[q,c+]_ := qc+ - c+Q = (15) 

respectively 

[c i ,c i ]_ = , [c+ c +]_ = . (16) 

These are exactly the raising and lowering operators as given explicitly in 
Brunei, Oerding and van Wijland (2000)1 . 



The dynamics is expressed by 

v / i=i 

where the master equation has to be inserted and evaluated to obtain the specific 
form of the operator L. The Liouville operator is after some calculation 

N I N A N 

L = ^(l-c i )pl^J ij c+c j \4 + J2(l-ct)ac i . (18) 
i=i \j=i J i=i 

The term (1 — a) guarantees the normalization of the master equation solu- 
tion and /3JijCjCjcf creates one infected at site i from a neighbour j which 
is itself not altered. Cj c j is simply the number operator on site j. Fur- 
theron, eta removes a particle from site i, again ensuring normalization with 
(1 — cf). This form Eq. (|1U|) and Eq. (|18|) is exactly the form given as well in 
[Grassberger and de la Torre (1979)] , there pp. 392-394, Appendix A), hence 
using the raising and lowering operators. 

Mean values for any quantity can now be given with a Felderhof projection 
state [Felderhof (1971)| . 

2.1 Perturbation analysis gives critical threshold 

The Liouville operator can be given in the form of a perurbation ansatz with 
an easily diagonalizable free operator Wo, acting only on single sites, and an 
interaction operator V contributing with strength A to the interaction (see for 
the 1 dimensional contact process e.g. [de Oliveira (2006) | ) 

L = W + A • V (19) 

with 

N 

W :=J2^i ( 2 °) 

8=1 

with Bi := (1 — cf) a an d without loss of generality a = 1, and with 

JV 

V :=J2Qi(ni-i + n i+1 ) (21) 



8=1 
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where Qi := n(l — c i) c t with Q = 2 the number of neighbours in two dimen- 
sions, and finally originating from the interaction term in one dimension 

N 

Yl J H C 1 C 3 = 4-1^-1 + 4+i c i+i ■ ( 22 ) 

J=l 

Identifying 

A := ■ Q (23) 

completes the perturbation ansatz Eq. (|19p for the Liouville operator given 
in Eq. (|18p . as given in de Qliveira (2006)] . With this ansatz the critical 



threshold and the critical exponents can be calculated very accurately via a 
scaling argument, e.g. for the time correlation function using the spectral gap 
de Qliveira (2006)] , and Pade approximation. 



3 Generalization to the SIRI epidemic model 

To generalize from the SIS epidemic model to more general cases like the SIRI 
we need to extend the number of dimensions of the basic vectors for the states 
(see e.g. Hinrichsen (2000~)] , p. 132), hence one more dimension for another 



particle type R besides I. For SIS we have 

|0> : = ( i ) = : l 5 > > I 1 ) : = ( I ) = : l J > ( 24 ) 
For the SIRI model hence we generalize to 



\S) := , \I):= \ 1 \ , \R) := . (25) 






Then creation and annihilation operators are given by 

a+\S) = \I) , a\I) = \S) (26) 
for infected, and for recovered we have 

b + \S) = \R) , b\R) = \S) (27) 
being achieved by the following matrix representation 

D \ 

(28) 
























1 


1 . -1 










V o 









V o 


1 


o / 



and 



2 
























1 ' -1 








° 


V o 





o J 




V i 





o / 



(29) 
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We can include the old two dimensional creation and annihilation operators, 
now extended to three dimensions, as 



6+ a , c := 1 = a + b (30) 





which can be expressed in terms of a + , a and o + , b. They give 

c + \I) = \R) , c\R) = \I) . (31) 

These operators a + , a, 6 + , b and c + , c are the ladder operators of the 8 Gell- 
Mann matrices X%, A 2 , ••• , As 

a+ = i(A 6 + fA 7 ) , a = ^(A 6 -iA 7 ) , (32) 

o+ = i(A 4 + U 5 ) , 6=i(A 4 -a 5 ) , (33) 

c + = i(Ai + »A 2 ) , c=i(A!-tA 2 ) , (34) 

with the Gell-Mann matrices (named after Murray Gell-Mann, a generaliza- 
tion to the Pauli matrices), well known in quatum chromo dynamics (QCD) as 
representations of the SU(3) group (special unitary group of dimension 3) 



A 2 = I i ) (35) 

(36) 

A 7 = I -i I (37) 

(38) 



and 





applied usually in QCD to the color states 

\red) = f I , \green) = ( 1 | . \hluc) - I I 






See e.g. [H iesmayr (2006)] and Koniorczyk and Janszky (2001)] . 
For the commutator rules we now have 

a a + = 1 - (a + a + 6+6) , (40) 

and 

b b + = 1 - (a+a + 6+6) , (41) 
see |Park and Park (2005)| . 
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3.1 The SIR! epidemic model 

We consider the following transitions between host classes for N individuals 
being either susceptible S, infected I by a disease or recovered R 



S + I 
I 

R + I 
R 



I + I 
R 

I + I 

S 



resulting in the master equation van Kampen (1992)] for variables Si, 1% and 
Ri £ {0,1}, i = 1,2, N, for N individuals eventually on a regular grid, with 
constraint Si + 1% + R% = 1. 

The first infection S + I — — > I + I occurs with infection rate j3, whereas 
after recovery with rate 7 the respective host becomes resistant up to a possible 

reinfection R + I — ► I + I with reinfection rate j3. Hence the recovered are 
only partially immunized. For further analysis of possible stationary states we 
include a transition from recovered to susceptibles a, which might be simply due 
to demographic effects (or very slow waning immunity for some diseases). We 
will later consider the limit of vanishing or very small a. In case of demography 
that would be in the order of inverse 70 years, whereas for the basic epidemic 
processes like first infection (5 we would expect inverse a few weeks. 

The master equation is explicitly given in the following form, as e.g. de- 
scribed in Stollenwerk, Martins and Pinto (2007)| 

— p (Si, Ii, R\, S2, h, R2, Rn, t) 



N 



N 



Y,P \J2 Ji ^i P{Si,h,Ru-A- Si,l- Ii,Ri...,R N ,t) 



i=l 



N 



+ ^7(1-/0 P (S 1 ,Ii,R 1 ,...,S i ,l-I i ,l-R i ...,R N ,t) 



(42) 



N 



N 



+ J2& \J2 3 vh p(Si,Ii,Ri,...,S i ,l-I i ,l-R i ...,R N ,t) 



1=1 

N 



+ £a(l-i? i ) p(Si,h,Ri,...,l - Si,Ii,l- Ri...,R N ,t) 



i=i 

N 

E 

i=i 



N 



V 



P Yl Ji i T i \Si + iIi + P\ Y. Ri + aRi 



J =1 



■p(...Si,Ii,Ri...) 
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J% j £ {0, 1} are the elements of the N x N adjacency matrix J, symmetric and 
with zero diagonal elements. The formulation of the master equation is given 
in analogy to the one used e.g. in Glauber (1963)| . 

3.2 Transitions in the SIRI model 

From the rules above we can obtain the different parts of the Liouville operator 
corresponding to the transitions in the SIRI model (derived from the master 
equation transition rates). 
Explicitly we have 

S + I^I + I , a +\S) = \I) (43) 

means creation of an infected from a susceptible in interaction with another 
infected, hence 

N / N \ 

Lp = P J> " Oi) £ j «j 4 ■ (44) 
i=i \j=i J 

Next we have 

R S , b\R) = \S) (45) 
means creation of a susceptible from a recovered, hence 

JV 

L Q = a^(l-6+)6, . (46) 
i=l 

Then we have 

I^UR , c+|/) = |i?) (47) 
means creation of a recovered from an infected, hence 

N 

i=i 

N 

= 7^(1 -0+^)6+0* (48) 
i=l 



Finally, we have 



R + I^Ul + I , C \R) = \I) (49) 



means creation of an infected from a recovered in interaction with another 
infected, hence 



N IN 

1 ij<-j<-j 



i=i \j=i 

N / N 

p - 6+ ai ) ^ Jijapjbfaj ] o+6j (50) 
i=l \j=i 



8 



With the commutation rules, especially Eq. (|4ip . we obtain 



c c 



,+ 



a + b b + a = a + (t — a + a — b + b)a = a + a 



(51) 



the number operator for the infected, hence 



H = 



N I N 




(52) 



Then the total Liouville operator is given by the sum of the individual 
Liouvilleans for the different transitions 



where the first two terms L 7 + L a are single site or free contributions and the 
last two Lp + Lp interaction contributions to the Liouville operator, suitable 
for a perturbation analysis analoguous to ansatz Eq. (|19p , 

4 Conclusions 

We have formulated the stochastic spatial reinfection model SIRI in terms of 
creation and annihilation operators in a Schrodinger equation like form, which 
opens the way for further analysis e.g. in terms of series expansions via a 
perturbation ansatz to calculate critical thresholds and critical exponents. The 
creation and annihilation operators turn out to be the ladder operators of the 
Gell-Mann matrices of the SU(3) group representation. 
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